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Abstract 

We perform traction experiments on viscous liquids highly confined between parallel plates, 
a geometry known as the probe-tack test in the adhesion community. Direct observation dur- 
ing the experiment coupled to force measurement shows the existence of several mechanisms 
for releasing the stress. Bubble nucleation and instantaneous growth had been observed in 
a previous work. Upon increasing further the traction velocity or the viscosity, the bubble 
growth is progressively delayed. At high velocities, cracks at the interface between the plate 
and the liquid appear before the bubbles have grown to their full size. Bubbles and cracks 
are thus observed concomitantly. At even higher velocities, cracks develop fully so early that 
the bubbles are not even visible. We present a theoretical model that describes these regimes, 
using a Maxwell fluid as a model for the actual fluid, a highly viscous silicon oil. We present 
the resulting phase diagramme for the different force peak regimes. The predictions are com- 
patible with the data. Our results show that in addition to cavitation, interfacial cracks are 
encountered in a probe-tack traction test with viscoelastic, liquid materials and not solely 
with viscoelastic solids like adhesives. 

1 Introduction 

Some materials display immediate stickiness, a property known as "tackiness" ^ E]- The de- 
formability of such materials enables them to achieve a good contact with all kinds of solid bodies, 
including those with surface roughness: Dahlquist's criterion jSlE), widely used since the mid- 
1960s, states that a solid material is sticky if its elastic shear modulus is lower than lO^Pa. In 
order to be usable as an adhesive, a deformable material should not flow on large time scales: it 
must be a viscoelastic solid. 

In a classical test called the probe-tack test in the adhesion community, a thin film of 
adhesive material is deposited on a planar, rigid surface. It is then tested with another planar, rigid 
surface, while the force is being recorded. The film is first compressed. After some "contact time" , 
it is subjected to traction. The traction force displays two characteristic features |S] before the 
separation is complete: a peak and a plateau. The reason for the force being relieved immediately 
after the peak can be traced back to two main mechanisms in usual adhesives. The most common 
one is cavitation, as evidenced by the first tests that included direct visualisation through the 
sample thickness Another, classical relief mechanism is the propagation of interfacial cracks. 
This occurs especially [5| when the elastic modulus of the material is high (low temperature or 
dense cross-linking). 

Whether the cracks appear at the sample edge ("external cracks" ^) or on multiple spots at 
the sample/indenter interface ("internal cracks") is mainly a question of sample aspect ratio [7||H]. 

^We use the names introduced by Crosby et al.T7\ 
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We here concentrate on how "internal" cracks and cavitation compete and interact. To address 
these questions, we continue the approach we used for cavitation [D] and study cavitation and 
crack phenomena in model material (viscoelastic liquids) both from an experimental and from 
a theoretical point of view. On such liquids (silicon oils), observed failure mechanisms include 
cavitation [T^ as well as "external" cracks from the sample edge m] in JKR geometry. 

In the present work, we first reexamine the usual cavitation and crack criteria (section . We 
then describe the protocols and materials used (sectional) . We present the experimental results and 
observations and offer a description of the underlying mechanisms (section 0J). We then construct 
a theoretical model to account for the triggering of the observed mechanisms (section ISJ and 
compare its predictions with the experimental results. We finally provide some discussion on the 
compared rheology of silicon oils and real adhesives in the context of our experiments (section (SJ. 

2 Cavitation and crack thresholds 

Cavitation and crack mechanisms are commonly encountered in adhesive films under traction. In 
the present section, we review and discuss the threshold stress needed to trigger them in the case 
of a purely elastic material. 

2.1 Cavitation and crack: an introduction 

2.1.1 Cavitation 

Cavitation in elastomeric materials under traction has been known since the experiments and 
calculations by Gent et al in the 1960s 221- The corresponding threshold for cavitation reflects 
the elastic resistance that the material opposes to the growth of inner, preexisting bubbles, and 
is commonly used in the context of adhesive materials [2] . For a full description of the cavitation 
process, one needs to consider other physical ingredients which affect the pressure required for 
bubble growth: the dilation of the bubble gas during bubble growth and the corresponding lower 
pressure, and the bubble surface tension which tends to make it shrink |13| . 

2.1.2 Crack 

Apart from cavitation, another mechanism is commonly encountered in adhesives films under 
traction: cracks often develop at the interface between the adhesive film and the indentcr. 

2.1.3 Method 

In the present section, we consider an elastic material that initially contains nuclei for both mech- 
anisms: bulk microbubbles (radius Rq) and microscopic cracks (size h) at the interface with the 
indenter. We take into account the possibility of propagation of interfacial cracks as well as all 
three ingredients involved in cavitation (elasticity, surface tension, gas pressure), and provide a 
very crude formulation for the corresponding stress threshold, restricting ourselves to scaling laws. 

2.2 Cavitation 

We now examine the physical ingredients that determine the cavitation threshold, discussed by 
Gent and collaborators |12l llHj . for an elastic material initially containing microbubbles. see Fig- 
ure ^ 

2.2.1 Elasticity 

In the regime where the cavitation threshold essentially reflects the elastic resistance of the material 
to bubble growth, the critical stress was calculated by Gent et al in the 1960s jj^l. In the case of 
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Figure 1: Left: bulk cavitation. Three main physical ingredients resist cavity growth and de- 
termine the (non-homogeneous) cavitation pressure threshold in an elastic medium: gas dilation 
(contribution on the order of the atmospheric pressure), cavity surface tension 7 and elasticity 
(modulus G) of the medium that surrounds the cavity. 

Figure 2: Right: interfacial crack. Main physical ingredients that determine crack propagation 
at the interface between a solid body and a deformable elastic material: applied stress a, elastic 
modulus G, local separation energy W for propagation (equation (T)), atmospheric pressure Patm- 
The dimensions h and 5 of the crack are also indicated. 



a neo-Hookean material, it is on the order of the (shear) elastic modulus: 

CGont — G (1) 

2.2.2 Dilation 

The growth of a microbubblc to millimetric size implies a strong dilation of the enclosed gas. 
The bubble growth rate is usually by far too fast for any gas diffusion from the bulk towards the 
growing bubble to develop significantly. As a result, the pressure in the growing bubble drops 
by an amount that is on the order of the atmospheric pressure. This contributes towards the 
cavitation threshold stress. When this term is dominant, the threshold is therefore: 

f dilation — Patm (2) 

2.2.3 Surface tension 

The surface tension at the bubble interface also contributes towards the cavitation threshold stress. 
When this term is dominant, the threshold therefore reflects the Laplace tensile stress exerted by 
the bubble interface on the elastic sample outside the bubble. It is proportional to surface tension 
and to the curvature of the bubble surface: 

f Laplace — "jT' (3) 

2.2.4 Cavitation threshold 

All three above ingredients enter the cavitation threshold for an elastic material initially containing 
microbubbles. A rough, simplified expression for the threshold is obtained as the sum (or the 
maximum) of all three values: 

<-^cavitation — Patm G -\- — — 

~ max|patm;G; (4) 



The value of these thresholds is reported on figure El 
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Figure 3: Gent- Tompkins diagramme for the stress threshold values for a bulk spherical cavity 
(log- log plot). The effective threshold is the largest of the three values given by elastic deformation 
(G), gas dilation (patm) and Laplace pressure {j/Rq). 

2.3 Interfacial crack propagation 

Let us now assume that disk-shaped cracks of size b are present at the interface between the 
(elastic) adhesive material and the (undeformable) indenter (see figure ISJ- We are interested in 
the value of the tensile stress that is required to induce the propagation of such internal cracks^. 

2.3.1 Thermodynamic work 

The thermodynamic energy Wq involved in opening such a crack includes the surface tension of 
the destroyed (adhesive/indenter) interface as well as those of both newly created (adhesive/air 
and indenter/air) interfaces: 

Wo = 7adli + 7ind - 7ind-adh (5) 

This Dupre energy is the simplest version of the work needed to separate both surfaces on the 
molecular scale. 

In practice, the energy needed locally to detach the adhesive from the solid substrate is larger 
than Wo- More elaborate estimations include local dissipation mechanisms such as the role of 
polymer molecules at or near the interface |14l I15| : 

W > Wo (6) 

In general, the interfacial energy cost corresponding to the crack of size b can thus be estimated 
as 

W b^ (7) 

2.3.2 Griffith's crack propagation criterion 

When a uniform, normal, tensile stress a is exerted onto the elastic material, the presence of the 
interfacial crack induces a slight reduction of the elastic energy since the crack essentially cannot 
transmit stress but is able to provide some extra volume to neighbouring regions. If the crack 
width b is increased, the interfacial cost (equation0) is enhanced while the elastic energy is further 
reduced. For a high enough value of the applied stress a, increasing the crack width b reduces the 
elastic energy to a greater extent than it increases the interfacial energy. As a result, the crack 
propagates under such a high applied tensile stress. 

This condition for crack propagation |1Q is known as Griffith's criterion for crack. Omitting 
numerical prefactors of order unity, it can be written as: 

y/GWjb (8) 

^We use the names introduced by Crosby et a/.|7|: internal cracks are located at the adhesive/indenter interface, 
while external cracks propagate from the edge of the adhesive / indenter contact region. 
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2.3.3 Crack and dilation 



GrifSth's approach was introduced in the context of liard, hardly deformablc materials. In such 
a context, the crack thickness 5 (see Figure is still very small at the onset of propagation, 
and the crack volume is thus always very small prior to propagation. For softer materials such 
as adhesives, the crack volume may increase sufficiently for the work done against atmospheric 
pressure to become predominant over the (Griffith) elastic and interfacial work. The propagation 
threshold is then on the order of Patm- As a result, the crack threshold can be reformulated as: 



f^surf — T^atm 



~ max-^Patm; y— ^ r (9) 

Expression where all numerical factors have been omitted, extends equation (jSJ) to softer 
materials or weaker interface strengths, which may be relevant in some cases for adhesives. The 
asymptotic regimes of this expression are presented schematically on Figure^ 
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Figure 4: Stress threshold values for a thin, interfacial crack. The effective threshold (expression^ 
is the largest of Griffith's value ^jGW/b and of atmospheric pressure Patm • 

We now compare the above expression for the interfacial crack threshold and the bulk cavitation 
threshold l@J. These thresholds will be central in our understanding of the experiments presented 
in sections 13 and 0] 

2.4 Competition between crack and cavitation 

Let us now determine how interfacial failure and bulk cavitation compete in the case of a purely 
elastic material. Since both the interfacial cracks of initial size h and the bulk cavities of initial 
size i?o (see Figures ^ and [2J) are subjected to the same applied stress cr, the failure with the lower 
threshold will trigger firsf^. 



2.4.1 General expression for the threshold 

Comparing expressions Q and lO, we can therefore approximate the global failure threshold as: 



cr ~ min [(Tbulk; Csurf] 



max<^Patm;G; 



i?0 



max i pati. 



WG 



(10) 



^Expression J^J is not always valid, however, as its derivation assumes that the shape of the crack remains 
disk-like. In other words, until propagation occurs, the crack thickness & must remain smaller than its width b. 
Also, the role of the trapped air IIYI . which partly relieves the pressure difference with the outside air, is not taken 
into account. A more elaborate discussion will be presented separately. 

''The initially triggered failure mechanism may not be predominant eventually, as discussed later in this article. 
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2.4.2 On the size of bulk microbubbles 



In view of Figure |31 tlie initial size Rq of bulk microbubbles is sometimes important to determine 
the cavitation threshold. 

In practice, in order to form small microbubbles (such that "f/Ro > Patm), one needs to 
incorporate small amounts of gas in the material during the formulation process. In order to 
achieve that, one needs to apply high enough stresses to overcome the corresponding Laplace 
pressure, i.e., stresses in excess of 10'"' Pa. 

In the following, for the sake of simplicity, we shall assume that this is not the case, i.e., that 
only larger bubbles are present (7/i?o < Patm)- Hence, the dilation contribution dominates over 
the surface tension contribution. 

As a result, the 7/i?o term in expression 110(1 can be left out, and the general expression for 
the threshold can be simplified as: 



max {patm ; G} ; max < patn 



WG 



(11) 



2.4.3 Large bulk microbubbles and strong interface 

Let us now consider an elastic sample with large bulk microbubbles (such that ^/Rq < Patm)- 
Equation is then especially interesting for a strong interface (W > fcpatm), as represented on 
Figure |S1 

• For very soft materials, both the cavitation threshold and the crack threshold are close 
to Patm- It is therefore impossible to determine simply which mechanism will occur (bulk 
cavitation or surface crack). It probably depends mainly on the local disorder in the material 
(G) or in the interface (W or b). 

• When the elastic modulus is increased (between letters A and B on Figure [Sj); the surface 
threshold (Griffith regime) becomes larger than the bulk threshold (at atmospheric pressure), 
and bulk cavitation is triggered first. 

• When the elastic modulus is further increased (between B and C), the bulk threshold remains 
lower than the crack threshold even though it now increases (Gent's regime), thus cavitation 
is still triggered first. 

• For large values of the elastic modulus (on the right-hand side of letter C), the cavitation 
threshold becomes larger than the crack threshold: surface cracks are then triggered first. 



2.5 Dahlquist criterion and Gent's cavitation threshold 

Dahlquist's criterion ^ for an elastic material to display adhesive properties states that its elastic 
modulus should be lower than around 10^ Pa. 

In view of the above discussion, given that this numerical value corresponds to atmospheric 
pressure, it appears that Dahlquist's criterion coincides with the crossover between two regimes 
for cavitation in an elastic material: Gent's elastic regime ((Tsouii = G) and the dilation regime 

(f^scuil — Patm). 

G~patm, (12) 

In practice, since the elastic modulus of many pressure-sensitive adhesives is lower than Patm, 
Gent's cavitation threshold may not be fully relevant for soft adhesives. Instead, we expect the 
dilation and Laplace cavitation thresholds as well Griffith's threshold for crack to be predominant 
in such materials. 

More precisely, if microbubbles are not too small (paragraph 12.4.2(1 and if the interface is 
strong, one expects the failure threshold for soft adhesives (G < Patm) to be always governed by 
dilation (atmospheric threshold, corresponding to the left-hand side of point B on Figure EJ- We 
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Figure 5: Stress threshold values as a function of the elastic modulus G (log- log plot) when bulk 
microbubbles are not too small {j/Ro > Patm) and when the interface is strong (W/(6patm) > !)• 
The overall stress threshold (thick line) results from the competition between bulk and surface 
thresholds. Hard materials (on the right-hand side of letter C) undergo interfacial separation 
(Griffith's criterion y/ GW/b). Moderately hard materials (between B and C) obey Gent's cavi- 
tation threshold G related to elasticity. Soft materials (on the left of B), which obey Dahlquist's 
criterion (i.e., G < Patm), all have an atmospheric stress threshold due to dilation. Moderately 
soft materials (between A and B) cavitate in the bulk, while very soft materials (G <C p^tm^/^' 
some distance on the left of letter A) exhibit either bulk cavitation or surface crack, depending on 
local material disorder. 



also expect the failure mechanism to be cavitation for moderately soft materials (between points 
A and B on Figure|3J), and either cavitation or crack for very soft materials (left-hand side of point 
A). 



3 Protocol, materials and experiments 
3.1 Apparatus 

The general geometry for a probe-tack test is the following: it consists in two horizontal and 
parallel plates whose separation h can be varied. One of the two plates is mounted onto a carriage 
(location t) via a load cell. The material is initially deposited onto the fixed plate. The moving 
plate is slowly approached, for instance until the material is confined into a film of prescribed 
thickness Hq. The material is then allowed to relax for a prescribed duration tc, known as the 
contact time. The carriage is eventually pulled at a constant nominal velocity V = i while the 
force F is being recorded. 

Two different traction machines are used for the experimental part of the present work. The 
first one is a commercial equipment (Z2.5/TN1S, Zwick Roell, Germany), the second one is a 
home-made prototype. The nominal separation velocity can be varied by about four orders of 
magnitude, from typically 1 /xm/s to lOmm/s. We usually mount load cells with a 100 N-capacity, 
but other transducers with a lesser capacity (for instance, 10 N) can also be used, if necessary. 
Both machines yield time t, force F and carriage location £ as digital data with an acquisition 
rate fixed at 50 Hz (commercial machine), or adjustable up to 1000 Hz (home-made prototype). 

A piece of polished, optical glass (BFI Optilas, France) is used as the fixed plate with the home- 
made machine. This allows to observe the bulk of the material during the traction experiment, 
and digital pictures (up to 1000 frames per second) may be recorded via a fast CCD camera 
(MotionScope lOOOS, Redlake,USA) with an optical field and a pixel resolution (up to 480 x 420, 
1-byte) depending on the chosen acquisition rate. 
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The fixed plate in the commercial machine is a square piece of anodized aluminum alloy. The 
probe is cither microscope slide glass or stainless steel (machine tool adjusted). 

3.2 Materials 

In this study, we have used a non- volatile silicon oil provided by Rhodia Silicones (France): Rhodor- 
sil gomme AS 522 (in short G20M), with a nominal viscosity of 20. 10"^ Pa s. 

Rheological curves have been determined at room temperature using a controUed-stress rheome- 
ter (AR2000, TA Instruments, USA), in a cone-plate geometry (diameter= 20mm and angle= 
4degrees). Two types of experiments have been conducted: oscillatory experiments (in the linear 
regime) and steady state flow measurements. Results are presented on figures HandlHl 

Figure El shows the linear viscoelasticity experiment. The sample displays an elastic behaviour 
(G" > G") at short times and a viscous behaviour (G" > G') at long times. The corresponding 
characteristic time can be defined, for instance, by the value 27r/Tci.ossing of uj for which G' and G" 
have equal values. 

The Cole-Cole diagramme presented on figure [7| shows that the rheology of the sample displays 
a Maxwell behaviour up to a; ~ 0.032i?z. We have extrapolated the Maxwell behaviour and 
determined a characteristic time TMaxweii, a plateau modulus Go and a viscosity {rj = G" /uj) (see 
lines on figure Eland on table 

Figure (HI displays the viscosity during steady state flow. The low-shear viscosity is Tystat. — 
20100Pa.s. This viscosity is similar to the viscosity measured by viscoelastic experiment (TyMaxwcU — 
20770Pa.s) and is consistent with nominal viscosity r^nominai = 20000Pa.s. 

We have also tested some non-linear aspect of the G20M rheology. At high shear rate, the 
sample undergoes bulk cracks and it becomes impossible to measure the viscosity in the steady 
state flow experiment. However, assuming that the steady-shear viscosity can be deduced from 
the oscillatory viscosity ("Cox-Merz's" rule ^H]), we can infer strong shear-thinning, see figure (HI 
(we could measure only about 10% shear-thinning directly in steady-shear). Such behaviour had 
been encountered with lower molecular weight silicon oils [S]. For oil G20M, the critical shear rate 
for the onset of shear-thinning, determined as the crossover between the asymptotic scaling laws 
for 77(7), is 7c — 1.55IIz. 

Table ^ summarizes all above results. Silicone oil G20M obviously does not have one single 
relaxation time. In sectional we shall nevertheless model it as a Maxwell fluid to account for the 
presence of both elastic (at high frequencies) and viscous (at low frequencies) behaviours. 

??M^^(Pas) rMw(5) Go(Pq) Tcrois(s) ic{s^^) 
G20M 20770 6.7 3100 0.94 1.55 



Table 1: Rheological properties of silicone oil G20M. 



3.3 Data processing 

A traction apparatus is not infinitely rigid. In particular, the force transducer has a finite com- 
pliance. The carriage location (. and sample thickness h do not differ by just a constant: the 
difference between them depends on force, F . Similarly, the sample thickening rate h differs from 
the nominal traction velocity V ^ I (this was first pointed out with a system in JKR geometry , 
then observed also in a flat geometry |10l EU| ) . 
We assume a linear machine compliance, 1/K: 

h{t) = ho + Vt-^ (13) 

where the carriage location during the traction on a sample of initial thickness Hq is written as 
£{t) ~ ho + Vt. For this expression to be valid, the force experienced by the material must be 
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Figure 6: Linear, dynamic viscoelastictic moduli G'{uj) and G"{uj) of G20M silicone oil. Open 
circles (o) represent the elastic modulus G", and close circles (•) the loss modulus G" . Both curves 
cross at Lu ~ ijJcrossing = 6.68 rad/s. The full lines correspond to the moduli of a Maxwell fluid 
with a viscosity of rjujw = 20770 Pa.s and an elastic modulus Go = 3100 Pa. 




Figure 7: Cole-Cole diagramme for silicone oil G20M. The sample behaves like a Maxwel fluid up 
to w ~ 0.032 rad/s. 
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Figure 8: Viscosity of silicone oil G20M. Closed circles (•) represent the steady shear viscosity. 
Open circles (o) show the complex dynamic viscosity obtained from oscillatory experiments. They 
obey Cox-Merz's rule over the entire range of accessible steady shear (frequencies up to 0.9Hz). 



fully relaxed to zero before traction starts: we systematically choose a long contact time tc- There 
only remains a small, static capillary contribution to the force. 

The compliance of each machine has been determined previously 0: A' = 4.5 10^ N/m for the 
commercial machine and K = 2.5 10^ N/m for the home-made prototype. 

4 Results and interpretations 

4.1 Results 

In this study, the traction velocity is varied in the whole available range, from a few /im/s to a few 
mm/ s and the thicknesses range from 100/im to We have observed two types of curves: 

• for low velocities or large thicknesses, the force decreases regularly after the initial force 
peak; 

• for large velocities or small thicknesses, the force presents a peak, a plateau and a subsequent 
force drop. 

Figureinidisplays some curves with a regular decrease and some curves with a plateau. Besides, 
we can observe that the plateau value increases with velocity. 

Visual observation (see figure [TT|l shows that the transition between both types of force curves 
corresponds to the transition observed recently between fingering and cavitation mechanisms using 
less viscous silicone oils [H]. We won't discuss this effect in the present article. 

At even higher traction velocities, (see figure [Tn|l . we observe that the plateau length decreases 
with increasing velocity, until it disappears at very high velocities. Besides, as will described in 
more detail below, there is a transition in the amount of material that remains attached to the 
indcnter when separation is complete (em cohesive versus arf/iesiuefailure) . 

This behaviour, which differs from the the previous ones, may indicate the existence of yet 
another failure mechanism. We will now discuss these observations in greater detail. 

4.2 Mechanism identification 

FigurelT^shows three pictures taken successively during the traction of a silicone oil G20M sample 
with initial thickness 200^m on the home-made prototype machine. It also schematically presents 
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Figure 9: Traction curves for silicone oil G20M at low velocities: V = O.OOljTim.s ^ (t), V 
OmSmm.s-^ (□), V ^ O.Olmm.s-^ (a), V = Ommm.s-^ (o), V = OMmm.s-^ (♦), V 
0.07mm. s^^ (>). 




Figure 10: Traction curves for silicone oil G20M at high velocities: V = 0.2mm. s ^ (t), V 
O.bmm.s-^ (□), V = O.lmm.s^^ (a), V = 1.5mm.s-^ (o), V = 2mm.s-^ (♦). 
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Figure 11: Photographs of sihconc oil G20M taken with the home-made prototype in the course 
of traction. With traction velocity V ~ Q.02mm/s (left), viscous fingering is observed and the 
force curve decreases smoothly after the peak (not shown). At ^ = 0.05mm/s (right), cavitation 
is observed, as well as weakly developed viscous fingering. The force curve displays a plateau (full 
diamond data points (♦) on Figure El). 

our interpretation of the mechanisms observed. On the first picture, the sample appears as the 
medium grey disk. Traction causes the sample to retract (second and third photographs). As it 
retracts, it leaves a thin film of silicone oil on the plate. This appears as the dark annular region 
with constant outer radius. A bulk cavity is visible on the second picture (a white circle has been 
drawn around it for clarity). This cavity is expanding (see third picture). This is the cavitation 
mechanism observed previously [SJ. 

Finally, yet another region can be seen on the third picture. It appears much brighter than 
the cavity. Besides, although its appearance and growth corresponds to a small increase in sample 
thickness, it has a large surface area. This suggests that it is thin. Moreover, we have noticed 
that after the experiment is complete, it corresponds to a place where the probe is free from any 
silicon oil. These observations lead us to believe that this region corresponds to an interfacial 
crack between the probe and the sample. 

In summary, this example shows two different types of cavities. The first one develops in the 
bulk: this is genuine cavitation. The second one grows at the interface and remains flat: this is a 
crack. 

4.3 Force curve interpretation: from cohesive failure to adhesive failure 
4.3.1 Origin of the plateau 

Let us explain why the force curve displays a plateau, whether cavitation or cracks appear. 

In a previous study [H], we have demonstrated that the existence of plateau in the case of 
bulk cavities is due to the difference between the very low pressure inside the bubbles and the 
atmospheric pressure outside the sample. The force drop after the plateau was interpreted as the 
penetration of air into the cavities. 

The explanation is similar in the case of cracks: cracks do not contain any significant amount 
of gas, and they are isolated from the outside air by the presence of a silicone oil seal. 

Let us now describe how the seal forms in practice. Figure 1131 shows a force curve with a 
few photos taken at a high traction velocity, and the main stages of the unsticking process are 
schematically shown on figure El 

A few small bubbles appear first f figure IHb and picture 1 on figure IT^ . Small cracks appear 
next ffigure IT^ and picture 2 on figure EJ- Almost instantly, cracks then grow and merge into 
a unique crack ffigure ITU and picture 3 on figure [T5|l . The crack stops before it has reached the 
sample edge, thus leaving an annular region free of any bubbles or cracks near the edge. This is 
the airtight seal. Since the stress is now relieved, the bubbles shrink back and no further cavitation 
can be triggered. The plateau is free from any further events (see pictures 4 and 5 on figure El- 
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Figure 12: Cavitation, crack and airtight seal. These three pictures have been taken successively 
with the home-made prototype during traction (silicone oil G20M sample with thickncss=200^m). 
The corresponding drawings expound our understanding of the mechanisms that take place in the 
sample. A bulk cavity is visible on the second picture (a white circle has been drawn around it 
for clarity) . This cavity is expanding (see third picture) . A whiter region can be seen on the third 
picture. It corresponds to a place where the probe is free from any silicon oil after the experiment 
is complete. These observations lead us to believe that this region corresponds to an interfacial 
crack between the probe and the sample. Two circles have been drawn on the third picture. The 
annular region between them is a free from any bubble or crack. It therefore isolates the inner 
region from the outside air and plays the role of an air-tight seal. 

Hence, the seal (see figure IT^ isolates the cavities and cracks (at quasi nil pressure) from the 
outside air (at atmospheric pressure). The airtight seal is thus essential for the existence of the 
force plateau. 

The force remains important (plateau) as long as the seal resists the pressure difference. When 
the seal eventually breaks, air comes into the cavity formed by the crack, and the force decreases 
abruptly (time (6) on figure EJ- The force drop corresponds to a pressure drop of roughly one 
atmosphere, as observed previously with cavitation P). 

4.3.2 Plateau length 

Figure [TUl shows that the force plateau is shortened as the traction velocity is increased. This can 
be understood qualitatively as follows. 

As the force drop is caused by air penetration ^1 , we expect the duration of the plateau to 
depend on the resistance of the airtight seal. Now, figure [T51 shows that the width of the airtight 
seal is reduced as the traction velocity is increased, which explains that the airtight seal then has 
a lower resistance. This reason for the decreased width of the airtight seal could probably be 
explained by a detailed theoretical analysis of the crack propagation in the present experimental 
conditions. 

4.3.3 Adhesive versus cohesive failure 

The traction velocity not only affects the plateau length but also the nature of the failure between 
the adhesive and the indcnter. Indeed, failure is observed to be adhesive in the region where 
the interfacial cracks have propagated, while it is observed to be cohesive in the airtight seal 
region. Hence, when traction velocity is increased, the annular region with cohesive failure becomes 
thinner. At very high velocities, interfacial cracks propagate right to the sample edge and no seal 
forms. Correspondingly, there is no plateau and the failure is purely adhesive (see figure EJ. 

The type of failure thus evolves from cohesive to partly adhesive and finally purely adhesive 
as the traction velocity is increased. 
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Figure 13: Traction curve for silicone oil G20M at a velocity of traction V = Imm/s with the 
home-made prototype. The photos were taken at six times indicated on the curve. Cavities first 
appear (photos 1 and 2). Interfacial cracks then appear and propagate (from photo 2 onwards). 




Figure 14: Major stages during separation. Some small cavities grow (a). Small interfacial cracks 
develop (b). The cracks quickly propagate further and merge, thus relaxing the tensile stress 
around the cavities, which therefore shrink back (c). 
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Figure 15: Sample aspect after complete separation for three different traction velocities (from 
left to right: V = 0, 13mm/s, V = 0,4mm/s and V = 0.7mm/s). The width of the airtight seal 
area (dark gray), which isolates the interfacial crack region from the outside air, can be seen to 
decrease with increasing traction velocity. 
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Figure 16: Force curve for a high traction velocity {V ~ 2mm/s). No force plateau is to be seen. 
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4.3.4 Conclusion of observations 



Upon increasing the traction velocity with such highly viscous G20M silicone oil, crack was ob- 
served beyond fingering and cavitation. We described the crack appearance, growth and merging. 
At the sample edge, we observed the presence of an airtight seal after cracks have merged. The 
airtight seal isolates the crack region from the outside air. It is thus responsible for the observed 
force plateau. The width of the seal is observed to decrease as the traction velocity is increased, 
until it docs not even form at very high velocities. One can infer that its resistance to the pressure 
difference between the crack and the outside air is weakened when its width is dereased. This is 
then consistent with the shorter force plateau. Only the crack region undergoes adhesive failure. 
Hence, the sample failure type evolves from cohesive to purely adhesive as the traction velocity is 
increased. 

5 Model and discussion 

As announced at the end of section we model the sample behaviour as that of a Maxwell fluid 
with characteristic time t. In the present section, we derive the expected behaviour for such a 
fluid in the context of our experiment. 

We first establish the evolution equation of the sample under traction fsection FS.lll . Therefrom, 
we then determine the main possible types of evolutions, independently of the failure mechanisms 
fsection l5.2|) . We then introduce these failure mechanisms and discuss qualitatively how they may 
be triggered and how they orient the evolution of the system fsection i5.3(l . It appears that in order 
to account for all observed phenomena, one must include the kinetics of cavitation fsection l5.4() . 
We then briefly discuss the triggering and propagation of intcrfacial cracks fscction i5.5|l . Knowing 
the influence of cavitation and crack, we then establish the phase diagramme of the system in 
terms of the experimental parameters fsection l5.6|) . We finally compare and discuss the theoretical 
predictions and the experimental measurements fsection l5.7f) . 

The experimental variables are listed in table |21 below, including in their adimcnsional version 
which we use in the remaining of this article. 



h 


H = h/ho 


sample thickness 


t 


T ^tV/ho 


time 


F 


Khc, 


force 


dh 
dt 


TT dH 1 dh 
^ ,\T V dt 


top plate velocity 


T 


■J- V T 

ha 


Maxwell fluid relaxation time 


a 


K hi 


average tensile stress 


2a 


Kh^ 


maximum tensile stress 
(at the center of the sample) 


f^scuil 


^ — CTseuil 


failure threshold 


^cav 


rr IT Patm a^Q 

J cav — X h„ 


cavitation time 


^growth 


rp l?,TTr]Val 
-t growth — Y 


bubble growth time 


^crack 


J crack — Kha\J b 


crack time 



Table 2: Adimcnsional variables 



5.1 Evolution equation 

We here present the essential ingredients that determine the sample evolution. A complete calcu- 
lation is to be found in Appendix fXl 



16 



with non-dimensional version: 



where 



If the material was purely elastic and homogeneous, with shear modulus G, the force would 
read: 

F='-^athlG^-^ (14) 
^=Cci^ (15) 

Parameter Cei is the ratio between the machine stiffness and the material elasticity (in our case, 
Col — 74 for the commercial traction apparatus, with sample diameter 2ao = 9.5nim and sample 
thickness ~ 200yum). 

Equation 1151 describes a disk of elastic material under traction. Derivating it with respect to 
time T yields, when H — \ <^1: 

^-C,i§^ (17) 

If the purely elastic material is replaced with a Maxwell fluid (Gr ~ 77), it is shown in Appendix fXl 
that T is simply replaced with T + T jT: 



•^+^ = ^01-^ (18) 



where T is the non-dimensional Maxwell time (see Table 12). 

Note that neglecting the elastic term T and keeping only the viscous term T jT^ in Equation H18|l 
yields the Stefan equation |23 used in our previous model [S]: 

^ = C^ (19) 

where 

G = Geir=--^ (20) 

is the ratio between the machine stiffness and the resistance of the sample once it has turned 
liquid. 

The force expressed by Equation (|18|l is transmitted by the machine, which behaves like a 
spring: 

T=l + T-H (21) 

where 1 -1- T is the motor position, and H that of the upper plate. By combining equation 12 II as 
well as its time derivative, with equation 1181 one obtains the evolution equation for a disk of a 
Maxwell fluid under traction: 

T {c^x^ + H -ij = T =l + T - H (22) 

This is the central equation in our model. In the absence of any cavitation or crack, it yields 
the evolution of the sample thickness, H{T), and that of the force, T{T), starting with initial 
condition H = \ at T = Q. Graphically, both quantities can be visualized very simply, as their 
sum is the uniformly varying motor position (see Figure [T7|l . 

5.2 Main types of evolution 

Starting with initial condition = 1 at T = 0, we now study the various types of system 
behaviours, in the absence of cavitation or crack. Only in section [5.31 will we study crack and 
cavitation. 
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Figure 17: Schematic representation of the system behaviour. The motor position 1 + T and 
the sample thickness H{T) are plotted as a function of time. The difference between the motor 
position and sample thickness values is the machine elongation; it is proportional to the force 
that is transmitted through the sample, as expressed by Equation H21|l . 



5.2.1 Viscous regimes 

Two types of behaviours have already been described in our previous work, dealing with purely 
viscous liquids [Hj. Example solutions are depicted on Figure \Th\ Beside numerical solutions to 
Equation H22|l . analytical approximate solutions can be obtained for various stages of the system's 
evolution in order to provide insight into the system's behaviour. The corresponding analytical 
expressions for the main variables arc provided in Tabled) 

For low values of C = Cc\T (dotted curve on Figure [T5|) . the sample remains still for a short 
period of time (stage V3); it then follows the motion of the motor very closely (stages VA and 
V6). 

By contrast, for large values of C (full curve on Figure IT^ . the sample remains still for a 
longer period of time while the machine elongates and the force rises (stage 1^3); only later does it 
suddenly flow and relieve the force (stage V5) while the upper plate catches up with the current 
motor position; it eventually follows the motion of the motor closely (stage VQ). 

Beside these viscous regimes, the elastic components of equation H22(l have several consequences, 
which we now discuss. 

5.2.2 Elastic behaviour at short times 

When JF = 1 — H is much larger than J-/T, the evolution equation (|22() reduces to its elastic 
version (equation EJ. As shown in Appendix fXI this occurs — unsurprisingly — at rather short 
times (T < T). 

When this elastic behaviour at short times is taken into account, the succession of stages 
depicted on Figure IThI becomes richer, see Figure [T^ 

The viscous stages, labeled 1^3 to V6, now follow clastic stage El; correspondingly, the regimes 
discussed above are labeled i?1346 (for C ^ 1) and i?1356 (for C ^ 1). 

Also, a third route has appeared, labeled i?246, which consists in elastic stage E2 followed by 
viscous stages and V6. Furthermore, equation can become invalid if large deformations 
are reached while the sample still behaves elastically, {i.e., prior to time T ~ T). Let us now 
consider both issues: the existence of two distinct clastic stages at short times, and the possible 
onset of large elastic deformations. 
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Figure 18: System evolution when the sample is essentially viscous. Top: sample thiekness as a 
function of time. Bottom: schematic succession of stages. When C ^ 1, the machine elongation 
1 + T — H remains much smaller than the sample thickness H at all times: stage ^3 (at short 
times) is soon followed by stages V4 and V6, as illustrated by the dotted curve, obtained from 
equation H22|l with C = 1/3 and Cci = 5.2. By contrast, when C ^ 1, stage ^3 extends over 
a longer period of time and the machine elongates much further; stages V5 and V6 then relieve 
the traction force, as illustrated by the full curve, obtained with C = 9 and Cd = 27. These two 
regimes, studied in our previous work |5] devoted to Newtonian fluids, correspond to the limit 
C^i/C ^ cx) in the present context of a Maxwell fluid. See Table 01 for approximate analytical 
expressions for the main variables during stages to V6. 
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Figure 19: Various possible sequences of stages (from El or E2 to 1^6). Regimes i?1356, i?1346 and 
R2A6 are the three evolution scenarii for the Maxwellian system in the absence of any cavitation or 
crack. The case where large deformations are reached while the sample behaves elastically would 
require further assumptions and is not addressed in the present work. 



5.2.3 Compared machine and sample compliance 

At very short times, when the sample thickness is still close to its initial value [H ~ 1), one can 
combine equations H15() and (|21|l and show that the evolution of the thickness and force is linear: 



H -1 ~ T 
T ~ T 



1 


1 


+ Cel 




Col 


1 


+ Ccl 



These equations reflect the fact that the deformation induced by the uniform motor motion (dis- 
placement T) is shared between the sample [H — 1) and the machine (F), according to the ratio 
Cei of their respective elastic compliances. As a result, elastic stage El or E2 arises at short times, 
depending on the value of Cc\- 

In stage El (with Cd 2> 1), the machine is more compliant than the sample: the sample 
deforms at a much lower velocity than the motor velocity {H ^ 1, see Tabled in Appendix IbI 
Correspondingly, in regime i?1346 of Figure 1201 the curve H{T) starts with a horizontal tangent. 

By contrast, in stage E2 (with Cei ^ 1), the sample is more compliant than the machine and 
deforms almost at the motor velocity {H ~ 1), and the curve H{T) in regime i?246 of Figure EUl 
starts with a slope almost equal to one. 



5.2.4 Large elastic deformations 

The elastic behaviour is well described by equations (j23|l and 1)24(1 . which are linear, when the 
shear deformation of the material in the gap is small. The shear deformation in the material is 
maximal at the edge and near the plates. It is equal to the product of the relative thickening 
{h — ho) /ho and of the aspect ratio a/h: 

7-(i/-l)J (25) 

Large deformations are reached when 7 becomes of order unity. Since a/h is much larger than 
unity, this occurs when H — 1 is still small. Hence, we can use the initial aspect ration, ao/ho in 
equation l(^ . 
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Figure 20: System behaviour depending on experimental parameters T and Coi (see equations El 
and ll8|l . Regimes i?1346 and i?1356, obtained when Cci/T ^ ao/ho and Cd ^ 1, result from the 
viscous behaviour of the sample coupled to the machine compliance, as described for a Newtonian 
fluid 9 . In the present situation of a Maxwell fluid, the succession of stages in the system 
behaviour is richer (see Figure ITiUl . In regime i?246, obtained when the machine is very rigid, the 
sample thickness follows the motor motion almost exactly. Finally, for large T (see equation I26|l . 
the sample reaches large deformations while still elastic. This regime is beyond the scope of the 
present work, as it would require additional assumptions concerning the mechanical properties of 
the material. 

Equation 1231 indicates that whenever 



large deformations are reached before the material flows {i.e., while T < T). 

Large elastic deformations cannot be addressed in the framework of linear elasticity. Treating 
them would require additionnal hypotheses on the material behaviour, which goes beyond the 
scope of the present article. The regime where such large elastic deformations occur is indicated 
on flgure EOl 

5.3 Triggering the failure mechanisms 

We have now determined the evolution of the system from equation H22(l . i.e., in the absence of 
cavitation or crack. The results are summarized in Figures [HH and HfH and in Table 01 

Let us now use the results of Section |21 to determine which regions of Figure 1201 correspond to 
cavitation or crack. 

5.3.1 Pressure as the triggering variable 

As discussed in Section |21 the relevant variable to determine when cavitation or crack should 
develop is the (tensile) pressure contribution due to traction. Since it is non-homogeneous in the 
sample, we take the highest value in the sample, which is in the center of the sample and equal to 
twice its average value. It is therefore equal to 2F/{il/h). In non-dimensional form, as indicated 
in Table m it is given by the product 2TH . Cavitation or crack is expected to develop when the 
value of this product exceeds the corresponding threshold determined in Section [21 



r > ^(1 + Cel), 



(26) 
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Figure 21: Sample thickness H{T) (left-hand side) and tensile stress J-H{T) (right-hand side) in 
regime i?1356 (see Figure EUl . with stages E\ (elastic), Vi (viscous), Vh (catching up) and VQ 
(weak force and constant velocity flow). 

The value of 2!FH{T) can be determined by solving differential equation (|22|l numerically. In 
order to determine the main regimes, however, it is sufficient to consider the expressions of !FH[T) 
during the various stages, which are given in Table H of Appendix IH) . 

5.3.2 Importance of the failure mechanism kinetics 

As a first approach, one might assume that whenever cavitation or crack is triggered, due to 
sufficient tensile stress, it relaxes the stress instantaneously. The experimental results presented 
on figure El as well as our earlier study [H], would seem to justify this assumption. 

However, our recent observations (see parts l4.2l and l4.3|l show that above some traction velocity, 
cracks appear after cavitation has started. This observation has two implications: 

• the cavitation threshold is lower than the crack threshold (indeed, cavitation appears first); 

• the stress relaxation induced by the cavity growth is not instantaneous. 

To understand these observations, we therefore need to take into account the kinetics of the cavity 
growth and determine its consequences on the stress evolution in the sample. 

5.3.3 Decorative versus effective cavitation 

Before we determine their growth rate, let us emphasize the fact that the cavity growth has two 
main consequences: 

• the cavity soon becomes visible (once it is around one micron in size); 

• the cavity later has a mechanical effect on the system (once its size has become comparable 
to the sample thickness, i.e. around one hundred microns). 

Since the required sizes for visibility and for mechanical effectiveness are very different, it 
may happen to be relevant to consider the period of time when the cavity is visible though not 
mechanically active. This stage is then called "decorative cavitation", as illustrated on figure 

5.3.4 Paths towards failure 

Depending on the growth kinetics, "decorative" cavitation may or may not play a significant role 
in the development of failure mechanisms. As a result, if we consider for instance viscous regime 
V?> as a starting point, various failure regimes can expected, as described on Figure 1^ 

If the crack threshold is lower than the cavitation threshold, only cracks are observed (regime 
m39). 
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Figure 22: Visibility and mechanical role of bulk cavities as they grow from their initial size Rq 
to their final, macroscopic size. As soon as their size exceeds about one micron, they become 
visible. They are not mechanically effective in relieving the tensile stress, however, until their size 
becomes comparable to the sample thickness. In the mean time, they can be adequately described 
as purely "decorative" . 

Otherwise, cavities appear and reach a visible dimension (decorative stage V7). The fate of 
the system then depends on the cavity growth rate: 

• If viscosity is low, cavities grow very quickly and cavitation can be considered instantaneous 
(regime i?138). 

• If viscosity is high, cavities remain decorative for a long time and effective cavitation is 
delayed (regime i?1378). 

• If viscosity is even higher, cracks may be triggered (stage C9) while cavities are already 
visible (regime i?1379). 

Let us now study in detail the kinetics of the cavity growth. 
5.4 Kinetics of cavitation 

Cavitation in an infinite and purely elastic medium was described by Gent and collaborators |12[ 
I13| . Recently, this approach was extended to the case of a finite sample [221123 to determine the 
final cavity size. 

In the present article, we essentially study cavities that appear in a viscous sample. Cavitation 
kinetics will thus be addressed only in the viscous stages described above: V3, VA, V5 and V6 
(see paragraph l5.2.1(l . As mentioned elsewhere jH), the kinetics of the microbubble growth (initial 
radius Rq) is governed by equation |24j : 



In this equation, p{t) is the pressure in the sample as determined in the absence of cavitation, 
and pcq{R) is the pressure at which a bubble of radius R is in mechanical equilibrium with its 
surroundings: 



The shape of function Pcq{R) yields a pressure threshold |5] (see figure [21 for an illustration) 
which corresponds to a tensile contribution slightly greater than Patm hi our case, as mentioned 
in paragraph 15.5.21 



R ^ Pcq{R)-p{t) 
R Ari 



(27) 




(28) 
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Figure 23: Various failure regimes expected to be triggered during viscous regime V3, taking into 
account the kinetics of cavity growth. If the cavitation threshold is lower than the crack threshold, 
cavities start to grow and soon become visible: they are decorative (stage V7). If viscosity is 
low (uppermost path), the cavities grow very quickly and become mechanically effective (stage 
C8); thus, the duration of the decorative stage is negligible and cavitation can be considered 
instantaneous (regime i?138). By contrast, if viscosity is high, cavities remain decorative for a 
long time and effective cavitation is delayed (regime i?1378). If viscosity is even higher, as the 
tensile stress continues to increase during the decorative stage V7, cracks may be triggered (stage 
C9). If cracks develop fast, they relieve the stress and hinder any further cavity growth. Hence, 
cavitation is observed for some time, but eventually cracks take over (regime i?1379). Finally, if 
the crack threshold is lower than the cavitation threshold, then the system evolves directly from 
V3 to C9. Only cracks are observed (regime i?139). 




Figure 24: Conditions of stability and growth of a bulk cavity, depending on the external pressure. 
The curve is a schematic representation of function peq(R) given by equation (|28|l . The evolution 
of the bubble size (illustrated by horizontal arrows) is determined by equation (|27|l . It implies 
that the first part of the curve (solid line) is a stable branch, while the second part (dashed line) 
is unstable. The corresponding critical point determines the quasistatic pressure threshold. 



5.4.1 Non-dimensional cavity grov^rth equations 

In non-dimensional form, equations H27|l and ()28|l read: 



d R - 

— = AR [2^F(T)-Seq(i?)] 
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(29) 



(30) 
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where R = R/Ro and: 



A 



3 1 



^ 



(31) 
(32) 



Here, Eoq is the non-dimensional form of Patm — Pcq, and C is the constant defined by equa- 
tion ||2n|) EI- Note the factor 2 in term 2TH{T) in equation It reflects the fact that the 
pressure (tensile) component due to the fluid flow is non-homogeneous in the sample and that in 
the center of the sample, it is equal to twice its average value. 

5.4.2 Cavity growth in the present experimental context 

As indicated in Tabled TH{T) is of order T in regimes El and Vi. 

The main trends of the cavity growth depend essentially on the initial slope B and on the 
maximum value Ec of function Y,c(^{R): 
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(33) 
(34) 



Note that B and Sc differ by a numerical factor which evolves in a limited range of values: 

3 < — < 3\/3 



(35) 



where the lower value, 3 corresponds to the limit 7 <C i?oPatm, and the greater value, 3'\/3, to the 
opposite limit, 7 » i?oPatm. 



5.4.3 Cavity growth parameter 

When solving Equation (|29|l for R{T), it appears that the dynamics of the cavity growth depends 
qualitatively on the value of parameter 



2 V^Tt/F 



27 



3/2 



3\/3 



27 



(36) 



This parameter is indeed relevant in our set of experiments. For the 77 = 10"^ Pa.s oil, with 
flo = 5 mm, Patm = 10^ Pa, K = 2 10^ N/m, taking a large traction velocity = 1 mm/s 
and assuming 7/i?o ^ Patm, one gets ScV^ — 1- Thus, as mentioned in our earlier work [H], 
the viscosity-delayed cavity growth appears even for an oil with viscosity 77 = 10"^ Pa.s at large 
traction velocities. A fortiori, the cavity growth parameter EcV^ takes small values with our 
more viscous oil [r] — 2 10* Pa.s). 



5.4.4 Instantaneous cavity growth 



At low traction velocities (ScV^ ^ 1), well before the cavitation threshold is reached (T = Sc/2), 
Equation ^ can be approximated as d.R/dr = AR [2TH-BR] by using = T and Ecq = BR. 
Hence, the cavity growth is mainly linear at short times, as illustrated on Figure !^ 



R~l 



2T 



-ABT 



AB^ 



2T 



(37) 
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Figure 25: Cavity evolution (time T, radius R) under low traction velocity (T,cVA » 1): instanta- 
neous cavity growth. In this regime, the cavity grows at first linearly (T < Ec/2), then much more 
rapidly (exponential growth). When the initial cavity is much too small to be visible (-Rvis S> 1), 
the observable growth (medium grey region) is thus very rapid (Tvis < T < Toft). By contrast, if 
the cavity is initially almost visible (^vis < 1 + ^c/B), then the duration of the observable growth 
is longer (hght grey), and a substantial part of the growth is essentially linear {T^^g <T < Ec/2). 



At later times (T > Ec/2), the cavity radius increases exponentially. 

As a result, the time Toff at which the cavities have a mechanical effect (see figure I22|) is 
essentially equal to Ec/2: 

As for the time Tvis at which cavities become visible (radius Rvis), it depends on how _Rvis 
compares with 1 + Ec/B, as illustrated on Figure!^ 

If the cavity is initially very small and thus needs to grow substantially before it becomes visible 
(^vis ^ 1 + Sc/-B), then Tvis is also on the order of Ec/2, and the cavity remains decorative (see 
figure 123) only very briefly: 

Tcff - Tvis = O (^-^^ « Y - - ^''^ ^^^^ 

By contrast, if the initial cavity is almost visible, i.e., 1 < R^is < 1 + Ec/i? (or already visible, 
with ^vis < l)i then it is possible to observe a slow, mainly linear cavity growth from T = Tvis to 
T ~ Ec/2, with Tvis ^ B{Ryis - l)/2 (or Tvis = 0, respectively). 



5.4.5 Delayed cavity growth 

At high traction velocities (Ec-s/A ^1), the cavity growth is exponential ^ and becomes substantial 
only well after the cavitation threshold (T = Ec/2) has been reached, as illustrated on Figure 1^ 

i?c.e+^^' (40) 
^This can be shown by considering inequality < Scq(-R) < Sc- Combined with equation 

dT _ 1 
dR ~ A R [2TH{T) - Scq(R)] ' 

it implies that ^^^^^ < ^< Aii(2T-Sc) ' '^^'^ hence that: 

logij < T(R) < Sc + logiJ 
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Figure 26: Cavity evolution (time T, radius R) under large traction velocity (Sc-\/A <C 1): delayed 
cavity growth. In this regime, the cavity growth is delayed by the fluid viscosity. The cavity radius 
increases essentially like e^^'^ . As a result, there is a substantial time lag ("decorative pause") 
between the time Tvis at which the cavity becomes visible and the time T^s at which it becomes 
mechanically effective (see figure 



In other words, the cavity growth is substantially delayed by the fluid viscosity, hence the term 
delayed cavitation. The time for visible and effective cavitations are then given by: 



logi?v 



A 



(41) 



Te« (42) 



i.e., typically: 
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(43) 



Toff ~ -= (44) 
VA 



In other words, in this regime of high traction velocity, the system marks a significant "decorative 
pause" and thus paves the way for a treaturous crack attack on cavity growth (see Figure 1^. 

5.4.6 Cavity growth during stage 

Among viscous stages V3 to V6, only stages V3 and V5 correspond to increasing tensile pressure 
(sec Table ]^ and are thus suitable for triggering cavitation. 

Yet, in the regime of delayed cavitation, once the threshold pressure has been reached (T > 
Ec/2), the system may evolve from stage V3 to stage VA (at time T ^ C during route i?1346) and 
still drive cavity growth. As the tensile pressure is constant during stage VA (with 2TH ~ 2C, see 
Table 21), the exponential cavity growth law is somewhat altered as compared to equation (|40|l : 

i?(T)^e+2^^^ = e^4^ (45) 
Hence, the times for cavities to become visible or mechanically effective become typically: 
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The conditions for this regime to arise arc the foUowing: 



C<~lMiogi?,ff <i (48) 

Ec/2 < C (49) 

Equation H48|l stipulates that the time at which cavitation is effective hes within stage VA, while 
Equation H49|l is the condition for the cavitation threshold to be reached prior to stage VA, i.e., 
during stage Vi. 

The second inequality in Equation (|48|l is always satisfied for thin samples. The other two 
inequalities can be satisfied only if: 



o ^2 

Sc < -^logi?cff (50) 
6 Oq 



^c<f^l0gi?eff (51) 



This is not the case in our series of experiments, since Cc is necessarily greater than 10^ Pa (see 
Section|21), while the right-hand side of if^ is on the order of 3 10^ Pa (with K ~ 4 10^ N/m, 
Co ~ 5 10-3 m and ~ 10"'' m). 



5.4.7 Elastic cavitation 

From the material point of view, the cavity growth implies a deformation mode (azimuthal stretch- 
ing around the cavity) that is distinct from the usual traction (shear in Poiseuille deformation 
towards the centre) . As a consequence, the resistance of the Maxwellian material to cavity growth 
depends on the growth rate. 

In the regimes described above, the cavity growth rate R/R is fastest when the cavity becomes 
effective. For cavitation developing from regime VZ, it is on the order of 2ATcs ~ 2^/A\ogR^. 
From regime VA, it is Za^/Ah"^. 

When R/R < l/r, the growth is liquid-like, as described above. When R/R > 1/r, however, 
the growth should depend on the elastic properties of the material. If the stress then exceeds the 
elastic modulus, i.e., Gent's threshold the cavity should expand to macroscopic (effective) 
size. By contrast, if the stress is lower than the elastic modulus, then the growth rate should 
stabilize at a value that allows viscous growth: R/R ~ 1/t. 

In the present situation, the elastic modulus (see Section I3.2|) is lower than the cavitation 
threshold (which is around atmospheric pressure). Hence, the cavity growth becomes elastic 
whenever the growth rate R/R exceeds l/r. 

5.4.8 When delayed is too late 

Cavitation from stage V5, called regime i?1358, implies that cavitation is instantaneous since stage 
V5 is very brief. In the present paragraph, we discuss whether viscously delayed cavity growth 
may hinder cavitation from this stage altogether. 

The duration of regime V5 is discussed in Appendix B.2 of Reference (Hj. It is on the order of: 

As for the maximum flow-induced tensile pressure component during stage V5, it was estimated 
as: 

•^^Imax - 4 - -g (53) 
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Using Equation (|29|l combined with both above equations, one obtains that during stage V5, a 
cavity can grow by a factor equal to at most 



e "0 



(54) 



where R^- (resp. Ry+) is the cavity radius immediately before (resp. after) stage V5. For the 
typical values of the sample dimensions ao — 5 10^'^ m and /iq ~ 10^^ m or ft-o — 5 10^^ m, one 
obtains: 



i?. 



R. 



{ho ~ 5 10"^ m) 



(55) 
(56) 



loading 



propagation 



Figure 27: Crack loading and propagation at the interface between a rigid body and a deformable 
material. The behaviour of a purely elastic material is well-known. At low tensile stresses (loading), 
the existing crack keeps its original dimension and the situation is quasistatic. Once the threshold 
stress has been reached, propagation occurs and the crack broadens very rapidly. For a viscoelastic 
liquid such as a Maxwell fluid, however, the (slow) loading stage (which is enabled by the hysteresis 
of the contact angle or by anchoring of the triple line) is followed by a quasistatic propagation 
(dewetting) which — let aside propagation velocity — is very similar to crack. In this stage, 
depending on the dynamics of the applied tensile stress, the dewetted region may widen at an 
increasing rate. Then, when the dewetting rate becomes high, the crack may behave elastically 
until propagation is complete. 



5.5 Crack triggering and propagation 

The question of the triggering and propagation of crack should be considered very carefully in the 
present context. 

Our material behaves roughly like a Maxwell fluid (see Section rO|l . and we have considered 
(see Figure I31|l that cracks are triggered while the macroscopic deformation in the material is 
viscous (T ^ T), at least in regimes i?1379 and i?139. 

The simple discussion of the crack threshold in paragraph l2.3l based on the assumption that 
the sample behaves elastically, must therefore be refined. This is particularly true for Griffith's 
criterion 

5.5.1 Elastic or viscous crack? 

In the present experimental situation, the applied stress increases linearly with time (stage El, 
E2 or 1/3). 

(i) At large loading rates, the material remains elastic until crack propagates, and Griffith's 
approach can be applied. 

(a) Conversely, at very low loading rates, the material behaves elastically only for a short 
period of time at early times. At all later times, it behaves as a liquid and may display dewetting. 
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Figure 28: Expected force peak time as a function of the traction velocity. Stage V3 (compliant 
machine and viscous sample) ends up with one of the following mechanisms: fingering, more rapid 
viscous flow {V5), instantaneous cavitation (C8, time Tcav), slow bubble growth {V7 between both 
dotted lines that represent times Tgrowth and 2rgrowth), or crack (C9 at time Tcrack)- Stage V7 
(slow bubble growth) can either continue up to the full bubble development (delayed cavitation), 
or be interrupted by crack propagation (the sample then displays both cavitation and crack). 



At moderate velocities, dewetting may resemble disk-like cavitation in the vicinity of the surface 
as described briefly in |25| . 

(Hi) At intermediate loading rates, as it acquires an increasing propagation rate, dewetting 
may progressively turn into elastic crack propagation. A detailed observation and analysis of such 
phenomena is given in Refs. [mi^l77| . 

5.5.2 Crack threshold in the present experimental situation 

It appears from the above discussion that in order to determine how interfacial cracks may be 
triggered and how they may propagate in the present context, a more elaborate discussion should 
be carried out and include the dissipation around the crack tip in a viscoelastic sample as it 
propagates [UlESIEZI- 

Such a detailed discussion goes beyond the scope of the present work. In the discussions below, 
we do not specify the expression of the threshold stress dcrack for crack. 

We are in a position, however, to provide some indications on the absolute magnitude of the 
effective crack threshold that should result from the considerations outlined above. Indeed, based 
on the observations of Section^and on the arguments of paragraph l5.3.2l it appears that interfacial 
cracks are triggered at a somewhat larger stress value than bulk cavities, and that the cavitation 
threshold is around atmospheric pressure. In other words, in terms of the (now obsolete) discussion 
on competing crack and cavitation in a purely elastic, solid material, the experiments reported 
here would correspond to a "moderately soft" material, i.e., located between points A and B on 
figure 

5.6 Summary: complete phase diagramme for fingering, cavitation and 
crack 

We are now in a position to predict the full system behaviour semi-quantitatively, and in particular 
the competition between cavitation and crack. 
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Regime A / Regime B 
Condition for Regime A 



Condition in terms of Cci and T 
Full, dimensional eondition 



i?246/i?1346 

Gel < 1 



i?1346/i?1356 
C < 1 



i?1356/i?1358 

C/8 < 



i?1358/i?138 
^/CJ2 < 



i?138/m378 



■logi?eff < Sc/2 



i?1378/i?1379 



ilogi?cfr < Scr/2 



i?1379/i?139 



i?139/i?19 
i?1379/i?19 

T < Ser/2 



i?1378/i?18 

r < 



i log i?cff 



i?138/i?18 
i?1358/i?18 

r < Se/2 



i?13* /L.E.D. 
T/Cc\ < ho/ao 



R246/L.E.D. 

T < ho/ao 



Ccl < 1 

^Gag < Khl 



CaT < 1 



CclT < 4Ee 



Patm 



32 h; 



-0 



^ V^V^A < aoCTcrack 



CelT^^gS2^./l0gi?v 

4Viogjj;^ie 



T < Scr/2 



-V^Vir < Op V Go-crack 



^rjVKho < alG 



2^J -K log R^a 



T < Sc/2 



T/Col < ho/ao 



\/t]V Kaoho < G 



T < ho/ao 
7]Vao < hloG 



Table 3: Equations that delineate the crossovers between the various regimes in phase dia- 
gramme EHl 
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Figure 29: Phase diagramme of the system behaviour in terms of non-dimensional parameters Coi 
and T (log-log plot). All three failure mechanisms (fingering, cavitation and crack) are included. 
The equations corresponding to all lines in the diagramme are to be found in Table |3 Varying the 
traction velocity as on Figure 1281 amounts to visiting a horizontal line on the present diagramme 
(since T oc V). 

This competition is summarized on figure which presents the time of the force peak as 
a function of the traction velocity in the case of fingering, cavitation or crack. More generally, 
the competition is illustrated as a phase diagramme on figure 1291 in terms of non-dimensionnal 
parameters Cd and T. Table |21 provides the equations for the crossovers between the various 
regimes in the phase diagramme 

We now review each regime very briefly and provide the predicted time for the force peak. 

5.6.1 Compliant sample and gentle flow: R246 

In this regime, which is achieved for instance for a rather thick sample (low Cd), the machine is 
more rigid than the sample. The sample therefore deforms and flows gently, almost exactly as 
prescribed by the motor. The sample eventually displays viscous fingering. 



5.6.2 Compliant machine and gentle flow: i?1346 

In this regime, as in most other ones below (Cd > 1), the machine is more compliant than the 
sample: it deforms more than the sample at early times. Here, the sample flows gently, almost as 
prescribed by the motor, and displays viscous fingering. This regime was described earlier [H] as 
"regime 1". 



'^pcak — 



I.e., 



tpoak ~ T ^ 1]/G 



(57) 



^pcak — C 




(58) 
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5.6.3 Compliant machine and sudden flow: i?1356 

In this regime, the sample resists traction for so long that it eventually flows in a very sudden 
manner (stage V5), after which it flows gently with the motor and displays viscous fingering. This 
regime was described in [0] as "regime 2" . 



Tpeakc^^^ *.e., Vak==^^/^;^ (59) 



5.6.4 Sudden flow and cavitation: i?1358 

In this regime, the very sudden flow induces a strong (tensile) stress peak (see Figure I21|l which 
triggers instantaneous cavitation. This regime was described in |^ as "regime 3" . 



Tpcakc^^y^ *.e., Vak^^^y^;^ (60) 

5.6.5 Instantaneous cavitation: i?138 

In this regime, the cavitation threshold is reached while the gentle sample flow is still unsignificant, 
and the cavity growth is so rapid that it can be considered instantaneous. 

Tpoak — S|;;/2 i.e., tpoak — 77 — (61) 

Z KV 

5.6.6 Delayed cavitation: i?1378 

In this regime, the cavity growth is delayed by viscous losses in the fluid. 



Tpcak - y ^ log i?off i-e., ipoak - \J '^Tl \og RcS \j (62) 

5.6.7 Decorative cavitation and crack: i?1379 

In this regime, the cavity growth is so much delayed that the tensile stress reaches the crack thresh- 
old and cracks propagate very rapidly. Meanwhile, however, the cavities have grown sufficiently 
to become visible, even though not enough to have any significant mechanical effect. 

Tpeak^Scr/2 I.e., W-^^crack (63) 

If the threshold stress (Tcrack does not depend on the traction velocity V , then the peak time ipcak 
is proportional to 1/V . 

5.6.8 Crack: i?139 

In this regime, cracks develop before the cavities could become visible. 

-^peak — ^cr/2 i.e., ^peak — ^^f^crack (64) 

5.6.9 Elastic cavitation: i?18 

In this regime, cavitation develops while the sample is still deforming as an elastic body rather 
than as a viscous material. This regime is not described here: the cavity growth implies large 
stresses and large local deformations. Additional assumptions on the material behaviour would 
be needed, and this goes beyond the scope of the present work. 
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Figure 30: Time of the force peak (upper graph) and appearance time of faihire mechanisms as 
observed with the home-made apparatus (lower graph) as a function of the traction velocity V 
(m/s). Circles indicate that viscous fingering was observed, and triangles correspond to cavitation. 
The sample used is the oil with viscosity rj = 1000 Pa.s. Power laws suggested by theory are 
indicated as guides for the eye. 



5.6.10 Elastic crack: R19 

In this regime, cracks propagate while the sample is still elastic. Again, this regime is not described 
in the present work. 



5.7 Test of the model against the experimental results 

Let us now compare our experimental results (see paragraph l4. II and our previous study 0) with 
the theoretical power law predictions concerning the time of the force peak (see figure |2H1 and 
paragraph I5.6|l . 

5.7.1 Experiments on the 1 GOG Pa.s oil 

Figure 1201 displays, among our previous results [^1, those obtained for the IGGG Pa.s oil. For 
comparison, two quantities have been plotted as a function of the traction velocity V (m/s) for 
these measurements carried out with the home-made apparatus: the time of the force peak (upper 
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Figure 31: Time of the force peak as a function of the traction velocity on the commercial (Zwick) 
machine The sample is the oil with 77 = 20000 Pa.s. Two series of experiments are presented 
(upper graph with open symbols and middle graph with filled symbols), which were conducted 
with two, in principle identical, steel indenters. Observation of the indenter after the separation 
is complete and interpretation of the traction curves indicate that fingering (circles), cavitation 
(triangles) or crack (squares) has occured. Power laws suggested by theory are indicated as guides 
for the eye. They arc reported on the lower graph for comparison. 
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graph) and the appearance time of failure mechanisms (lower graph) . Circles indicate that viscous 
fingering was observed, and triangles correspond to cavitation. The power laws suggested by 
theory for regimes i?1356 and i?1358 (sudden flow, slope — f/2), i?138 (instantaneous cavitation, 
slope —1) and i?1378 (delayed cavitation, slope —f/2), are compatible with the data, although 
not fully convinving, on the upper graph. They arc somewhat more convincing on the lower graph 
(observed appearance time for fingering or cavitation). 

The final —f/2 slope confirms the existence of delayed cavitation and refines the interpretations 
given in reference [Hj. 

5.7.2 Experiments on the 20 000 Pa.s oil 

The results obtained on the ?; = 20000 Pa.s sample are presented on Figure 1^ The time of the 
force peak measured on the commercial (Zwick) machine is plotted as a function of the traction 
velocity. Two series of experiments are presented. They were carried out with two, in principle 
identical, steel upper plates, which may differ slightly in — for instance — surface roughness. 
The symbol shapes indicate the failure mechanism that can be deduced from the observation of 
the indenter after the separation is complete and from the shape of the traction curve: fingering 
(circles), cavitation (triangles) or crack (squares). Crack corresponds to the absence of material 
on the indenter (adhesive failure). Cavitation corresponds to craters in the material that remains 
on the lower plate, or to the presence of a shouldering shape on the traction curve soon after the 
force peak. Fingering is revealed by visual observation of the material that remains on the plates. 

Power laws suggested by theory for regimes i?f358 (sudden flow and cavitation, slope —f/2), 
i?f38 (instantaneous cavitation, slope — f), i?f378 (delayed cavitation, slope —f/2) and i?f379 or 
i?139 (crack, slope — f) are indicated as guides for the eye. They fit the data rather convincingly. 

The straight lines for each indenter (upper and middle graph of Figure I3f |l are reported on 
the lower graph for comparison. It appears that at moderate traction velocities, when fingering 
(circles) or cavitation (triangles) occurs, the force peak time does not seem to depend on the 
indenter surface. Indeed, all data follow quite accurately the same power laws with exponents 
— f/2 and — f . Conversely, at higher traction velocities, the force peak times from each experiment 
are rather well described by a power law with the same exponent — f but with a different prefactor 
for each indenter. 

This behaviour (same prefactor) was expected for viscous fingering, which is a bulk phe- 
nomenon. As for cavitation, the fact that it does not strongly depend on the properties of the 
interface between the sample and the indenter seems to indicate that it may nucleate in the bulk 
(although we cannot draw a definite conclusion on this matter with only two indenters tested). 

Unsurprinsingly, as an interfacial phenomenon, crack is readily affected by the indenter surface: 
the threshold for crack propagation is observed to be somewhat lower for the second indenter than 
for the first indenter. 

In the case of the second indenter, the actual prefactors of the cavitation and crack regimes are 
very close to one another. As a result, the regime of delayed cavitation is not strikingly obvious. 
Only with the first indenter is it somewhat visible. 

5.7.3 On the onset of cavitation 

The reader may have noticed from Figure 1^ that all data points on the low- velocity —1/2 slope 
correspond to fingering: cavitation could not be infered from the shape of the force curve in this 
regime; and direct visual observation was not possible for these series of experiments (conducted 
on the commercial machine). 

In order to explain that no cavitation is present in the regime with a — f/2 slope, let us recall 
the discussion on delayed cavitation arising during stage V5 (see paragraph 15. 4.8(1 . 

In Reference for a 5 10~^m sample, the transition between i?1358 and i?138 is observed for 
C ~ 70 and the transition between i?1356 and i?1358 is observed for C ~ 1. The corresponding 
exponents in equation (|^ are 3.8 and 0.63, respectively. 



36 



On Figure 1^ the onset of regime i?138 (transition between slopes —1/2 and —1) is observed 
for V ~ 1.2 10~'^m/s, i.e., for C ~ 1.6 10"'^. As a result, the exponent in equation H54|) is 0.33. 

Hence, the cavity growth during stage V5 in the experiments reported here is expected to be 
less pronounced than in the experiments reported in Ref. 1^. This may be partly explain why 
cavitation is not observed in the regime where stage V5 is present (with slope —1/2). 

5.7.4 On the orders of magnitude 

As mentioned above, the expressions given in DaragraDh l5 . 61 account for the dependence of the peak 
time ipeak on velocity quite well (see Figure l?T|l . They do not, however, account for the correct 
orders of magnitude when the material parameters are taken as G = 3 10'^ Pa and ry = 20000 Pa.s. 

We believe that this discrepancy has its origin in the strongly non-Maxwellian character of the 
sample rheology (see Figures El and [Tjl . 

6 Conclusion and perspectives 

We conducted probe-tack experiments on highly viscous silicon oils. Beyond viscous fingering and 
cavitation reported in a previous work we observed delayed cavitation and interfacial fracture. 

We constructed a theoretical model of how a Maxwell fluid should behave in such a probe-tack 
experiment, including considerations on crack thresholds and on cavitation thresholds and growth 
kinetics. Meanwhile, we showed that atmospheric pressure contributes to the traction force both 
in the case of cavitation [51 ^1 and (in the present work) in the case of crack. 

Although the rheology of the silicon oils we used departs from that of a Maxwell fluid substan- 
tially, we were able to give a possible explanation for the existence of the various regimes observed 
experimentally: cavitation alone, delayed cavitation, cavitation followed by crack, and pure crack. 

Let us finally discuss two points. 

1. Why did our approach work at all? Why did observe phenomena not unrelated with those 
observed in true adhesive materials? 

2. What further rheological features should one include to mimic adhesive materials more 
closely? 

6.1 Why did this approach work? 

The above reported phenomena are very similar to those observed in adhesives. This may appear 
surprising, as the rheology of the systems we used (silicon oils) notably differs from that of adhesive 
materials In particular, silicon oils are viscoelastic liquids while adhesive materials are viscoelastic 
solids. 

In fact, this can be understood very simply by considering the possible rheological proper- 
ties |2H1 of a soft material (figure |221l- The distinction between solid and liquid appears only at 
long time scales: either the material develops a permanent resistance to flow (and it is a solid) or 
it eventually flows (it is then a liquid). At shorter time scales, only elastic and viscous characters 
are relevant. ^ 

Now, our theoretical predictions deal with two different stages in the course of traction: 

1. triggering the failure mechanisms; this happens while the material is still weakly deformed 
and still has not had time to display its solid or liquid character; 

2. the force curve after cavitation has occurred; the material has then been strongly deformed 
to allow for the cavity growth; the viscous description we gave of this deformation (where 
plastic may have been more appropriate) is qualitatively valid. 



^Indeed, beside the usual elastic solid and viscous liquid, soft materials include elastic liquids (generically 
represented by the Maxwell model) and viscous solids (generically represented by the Voigt or Kelvin model). 
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Figure 32: Rheology of soft materials under weak stress (schematic representation). On rather 
short time scales, depending on its molecular architecture (branching) and on frequency, the 
material may behave in a rather elastic or in a rather viscous manner. On long time scales, 
either it develops a permanent resistance to flow (and it is a solid), or it eventually flows (and 
it is a liquid). For some solid materials, a higher stress may trigger the flow (this is plasticity). 
Depending on the order of magnitude of the corresponding stress threshold, one usually refers to 
such a material either as a yield stress fluid or as a plastic solid material. 

6.2 Extending the material rheology 

The choice made in the present work to study a material with a Maxwellian behaviour (or more 
generally, any viscoelastic liquid material whose behaviour at large stresses is not specified) suffers 
some limitations. 

• Some (large traction velocity) regimes of the macroscopic sample deformation are not acces- 
sible (see region "Terra incognita" on Figure [501). 

• The large deformations around cavities can be adequately described in the slow, viscous 
regime but not in the faster, elastic regime. 

In order to address these questions, the material should either have a narrow elastic regime at 
weak stress, followed by a plastic behaviour, or be able to sustain very large elastic deformations 
before it yields (or hardens). 

True adhesive materials are often physically crosslinked and thus may display a plastic be- 
haviour, at least at moderate stresses. (At higher stresses and deformation rates, broken physical 
crosslinks may not have time to reconnect and the material may become thinner and eventually 
break). Thus, extending the material rheology to viscoelastic plastic solids and to large elastic de- 
formations will be important to capture more extensively the behaviour of true adhesive materials, 
essentially during three stages: 

• at the early stages of traction under large traction velocities, when large shear stresses 
develop within the confined sample; 

• once cavitation has been triggered, for the cavity development and growth kinetics; 

• once cavitation has fully developed, when the cavity walls experience continued stretching 
and may induce interfacial cracks 
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A Confined Maxwell fluid under traction 

In the present Appendix, we first compute the flow of a Maxwell fluid that is confined (thickness 
h, initial value ho) between two disks (radius a, initial value oq) when the disks are separating at 
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velocity h. We then apply the result to the probe-tack situation where the motor is pulling on the 
plates at constant velocity V via the force sensor which behaves like a spring. 

For later convenience, let us define the (constant) volume fl = ttoq Hq = ira^ h of the sample. 
Let also r be the distance from the axis of symmetry, and z be the altitude perpendicularly to the 
plates, with z = at mid-height z = ±h/2 at the plates. 

A.l Maxwell fluid in conflned geometry 

When the disks are pulled apart, the fluid is mainly sheared due to its strong confinement {h <^ a). 
Provided some approximations are valid (discussed in paragraph lA.SI below). the local constitutive 
equation for a Maxwell fluid can be written as 

where cr is the shear stress, e is the shear strain rate, G is the material shear modulus at high 
frequencies, and rj is the viscosity. 

A. 2 Velocity field 

As mentioned above, conflnement implies mainly radial velocity, in other words, the lubrication 
approximation is valid (except in the vicinity of the sample edges because of recirculation and in 
the very center where the radial velocity becomes smaller than the velocity along the vertical axis 

The magnitude of the radial velocity, averaged over the sample thickness, is fixed by volume 
conservation: 

+ h/2 

r ■ 

v{r, z) dz = v{r) = — h (66) 
■ h/2 2/i 

The velocity profile along direction z reflects the balance between shear stress and pressure gradient 
via the constitutive equation of the fluid. Classically, for a Newtonian fluid, the profile is parabolic. 
In the present case, since equation (I65|) is linear and since convective effects are negligible (see 
DaragraDh lA.5l below) . the profile is still parabolic: 



v{r,z) = v{r) ^ 



A. 3 Pressure field and total force 



ih/2y 



(67) 



The local stress balance drP — dz<7 (where a is the rz component of the stress) implies that the 
pressure gradient is related to the shear stress on the plates: 

^Iplatc = ^drP (68) 

Taking the time derivative: 

^Iple^te = ^drP + ^drP (69) 

Also, from equation (|^ . the shear rate at the plate is given by 

^Iplatc = ^ C^O) 



Combining equations ((HSJ, (EH), (EHJ and l(7n|) : 



3/i 
— r 

/l2 



h_ h_ 



drP{r) + ^drP{r) (71) 
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Taking p{a) = Patm and integrating from r to a, wc obtain an equation for the pressure field: 



3^1. 2 2n 



h h 



batm - p{r)] + - Pir)] 



(72) 



where p(a) = p{,'i^)\r=a ~ ^rPlr=a "^^^ ^'^ neglected according to assumption H79() . Integrating 
over the disk surface area and using F = [patm ^ pC^)] Sttj' dr, we obtain the equation for the 
force: 



3h 1 TT 

7^22 2 



2G ^ 2^ 



h ■ 

FH F 

2G 



Using Q = Tra^ h and -q = Gt: 



(73) 



2iT T \ h 



(74) 



A. 4 Coupling with the machine and evolution equation 

Using the spring equation 

F = K{ho + Vt^h) (75) 
and the adimensional variables (Table EJ, the differential equation (|74p can be written as: 



(76) 



The above equation, which is identical to Equation H22() . describes the behaviour of a Maxwell-like 
system in a probe-tack geometry. We now discuss its validity. 

A. 5 Validity of the local equation 

Equation H65() . which we used to derive the force response of the sample fEauation l74|l . involves 
the sole shear component arz of the stress. We shall now discuss whether or not it is valid to use 
this simple, scalar equation in the present context. 

The relevant tensorial equation for a Maxwell fluid such as a polymer melt is the uppcr- 
convected Maxwell equation: 



a'' + {v.\/)a'' - {VvY .a"" - a'Wv = Gi - 
where a'^ is the deviatoric {i.e., traceless) part of the stress: 

a-' = a- ^tr(a) 



(77) 



(78) 



The second term in Equation H65() is the usual gradient term in transport derivatives. The third 
and fourth terms, which involve the velocity gradient, and are one [upper- convected) form of the 
convective terms that are relevant when transporting a tensorial quantity that is linked to the 
underlying material medium. 

The use of the simpler equation instead of the full equation l(77)l implies that both following 
conditions be satisfied: 



(79) 
(80) 
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The components of these two tensor equations can be expressed as: 



Vrdra'^„ < af^ (81) 

Vr drCr^^ < (82) 
Vrdrat « (83) 



2cr^^ drVr + 2af^ d^Vr < cr,t (84) 
a^,. drv, + ai, d,vr < cr,t (85) 
2crf^ + 2crf^ a,.^^ < crf^ (86) 

Because the entire calculation carried out here is based on the lubrication approximation, the 
normal stresses <7„ and a^z cannot be distinguished from the hydrostatic pressure. As a result, 
the deviatoric normal stresses are zero: 

at, = 2at ^ (87) 

As a result, among the above six conditions fcouations 18 lKlSS|l . only Equation (|82|l provides a 
useable constraint: ^ 

a- 1 d fh^/^F\ 

/^•-•T^«3T (88) 



i.e., 

h F H T 

-h«F H«J 
From the material point of view, the stress must be weak enough for the recoverable deformation 
to be small: 

rr h^/'^F 

-^«1 I.e., (90) 



TH'''^ « Cei^ (91) 

Note that at short times, when the sample is elastic, weak stress implies small deformations, 
defined by equation (|25|) . In particular, it implies ~ 1. Condition (|91|l then reduces to: 

T « Cei^ (92) 
ao 

Applying this criterion to equation H24|l. one recovers condition (j26|l . 



B System evolution: stages and crossovers 

Here are two tables that summarize the results of the discussion in section 15.21 concerning all 
stages that can be encountered during a probe-tack experiment on a Maxwell fluid: 

• Table 0] provides the values of the main variables in all stages E\, E2, V3, VA, Vh and V(y. 

• Table m indicates the equations for the various crossovers. 

Table 0] describes only stages where the macroscopic sample deformation is involved. It does not 
describe other stages iVl, C8, C9), where the sample deforms locally around cavities or cracks. 
It provides the tensile stress TH , however, which is the relevant variable for triggering cavitation 
or crack (see paragra.ph l5.3.1|l . 

'^It is not excluded that other conditions provide stringent constraints when expressed beyond the framework of 
the lubrication approximation, but such a detailed hydrodynamic study is beyond the scope of the present work). 
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Stage name 


vaiiQiLy 


Main variable values 


El 


Ccl » 1 

T< r 

T < Cci ha/ao 




E2 


Cel< 1 

T<T 

T < ho/ao 


jc- ~ r Col 


V3 


Ccl » 1 

T> r 


J^ = T 

~ (1 - 2r2/c)-i/4 ^ 1 + ^ 1 

T H = T 




c< r< 1 


iJ = 1 + T - C ~ 1 + T ~ 1 




Cel » 1 


/cT 4\/2 

Y 2 ^ ^ p3/2 






2<H < 




C = CciT > 1 




V6 




J- H 



Tabic 4: Values of the main system variables dm-ing stages El, E2, V3, Vl, V5 and V6. 



El ^ 


E2 


Ccl - 


t 1 


machine and sample 










compliance competition 


El - 


V3 


T~ 


r 


Maxwell transition 


E2 - 


Vl 








V3 - 


V5 


H -1 


~ 1 


flow acceleration, force peak 


V3 - 


VI 


H ~ 


1 


flow stabilization, force peak 


V4- 


V6 


H-1 


~ 1 


sample starts deconfining 


V5- 


V6 


H - 1 


~ T 


fast to slow flow transition 



Table 5: Stage transition criteria and interpretation. 
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